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1. Introduction 


The SKA will observe the Epoch of Reionization (EoR) through the redshifted 21 cm emission 
of hydrogen. The intensity of the emitted signal depends on several local quantities but during most 
of the EoR the most important one is obviously the ionization state of the gas: only neutral atoms 
in the ground state can emit 21 cm radiations. Fully ionized regions will show up as null-signal 
areas in the tomographic data cube. While the basic physics of hydrogen ionization is simple, 
it takes place in a strongly inhomogeneous medium and thus involves vastly different scales and 
characteristic times, resulting in a complex geometry and evolution of the ionization field. 

To discuss the important aspects of the reionization process we will first consider that the gas 
consists of pure hydrogen. The effects of the presence of helium and metals will be discussed as 
needed. If x, is the local ionization fraction of hydrogen, it obeys the simple equation: 


dx f r°° <j I 

—^- = y(T)(l—Xi)n e — a(T)xin e + r(l—Xi) with T= dD. v v dv (1.1) 

dt J Jy Q hv 

where n e is the number density of electrons, T is the temperature, y is the collisional ionization 
coefficient, a is the recombination coefficient, T is the photoionization rate, v is the frequency, 
7 V is the local specific intensity, a v is the ionization cross-section, and Q. the solid angle. Tightly 
coupled to the ionization equation is the energy equation: 


dE 

dt 


Jf- A 


( 1 . 2 ) 


where E is the internal energy of the gas, Jff is the heating rate and A is the cooling rate. Both 
and A include contributions from radiative processes (photo-heating, collisional ionization cooling, 
recombination cooling, collisional excitation cooling, Bremsstrahlung cooling and Compton cool¬ 
ing - heating) and from dynamical processes (adiabatic cooling or heating, heating from viscous 
dissipation). References for fitting functions of the various rates, that depend on temperature and 
ionization state and for the ionization cross-section that depends on the frequency can be found for 
example in Iliev et al. (2006). 


2. Insights into the ionization process 

Along with the equations above, the process of reionization is governed by the radiative trans¬ 
fer equation (Rybicki & Lightman 1986). In the general case, these equations are furthermore 
coupled to the hydrodynamics equation of the self-gravitating gas through the energy equation and 
the equation of state of the gas. However ionization by UV photons will only heat up the gas to 
~ 20000 K, and the amount of X-rays produced during the EoR will not change this picture in 
the bulk of the IGM. The associated increase in the pressure of the gas will affect noticeably only 
structures with a virial temperature of the same order as the heated gas temperature or lower. If 
the radial density profile of real halos actually follow the numerical fit in Navarro et al. (1997) 
(the so-called NFW profile), structures with a virial temperature less than ~ 20000 K have a to¬ 
tal masse < 1O 8 M 0 (both dark matter and baryons included). The gas in these structures may be 
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photo-evaporated by the arrival of an ionization front or will never fall into the potential well of 
halos in ionized regions. On larger scales though, the coupling to the dynamics may safely be 
ignored as gravitation dominates over hydrodynamic pressure. 

Evaluating a few physical quantities in typical EoR environments will help us get a first idea 
on how reionization proceeds. 

2.1 Mean free path of ionizing photons 

In a homogeneous medium the mean free path of ionizing photons is defined by / = 

where a v ~ 6.3 10~ 18 ^ j end 2 is the ionization cross-section (see Verner et al. (1996) for a 
more accurate formula), Eq is the energy of the ionization threshold and //// the number density 
of absorbers, here neutral hydrogen atoms. In a fully neutral IGM at the average density of the 

universe we find Z ~ 2. (-pp)~ (p-) comoving kpc. In higher density regions the value will be 
even smaller. Consequently, a reionization fueled by stars that produce ionizing photons at energies 
not much higher than the ionization threshold will proceed in the form of fully ionized bubbles 
growing around the sources, separated from the neutral region by a sharp ionization front. On 
the other hand, the small amount of X-ray photons emitted above 1 keV have a mean free path in 
the 1 comoving Gpc range and will provide a weak (a few percent at most) and mostly uniform 
contribution to the ionization field (note that the associated heating has important consequences 
for the 21 cm signal). Photons with energies of a few keV have only a small probability of being 
absorbed in the IGM before reionization is complete. 

2.2 Recombination time 

Another important quantity is the typical recombination time. Or, in other words, how many 
times does an hydrogen atom in the IGM need to be ionized by the end of the EoR. The recom¬ 
bination time is defined as f rec = a( j jnil ■ At the average density of the universe and for a gas at 

T = 10 4 K, we find, t KC ~ 240 (j^) ^Myr. Since the bulk of the IGM is likely reionized in the 
6 < z < 10 range, most atoms need to be ionized only once. At the same time, a moderate over¬ 
density of ^ ~ 20 will have a much shorter recombination time of ~ 10 Myr. Such structures 
(filaments, minihalos) will act as photon sinks until the ionizing flux is large enough to photo- 
evaporate them (see, e. g. Sobacchi & Mesinger 2014). The recombination time in the interstellar 
medium will be even shorter and, most likely, only the dumpiness of the medium or holes blown 
by supernovae explosion or AGN feedback will allow ionizing photons to escape into the IGM (see 
next section). 

2.3 Ideal Stromgren sphere model 

Assuming sharp ionization fronts, neglecting recombinations in ionized regions and neglecting 
the coupling to the hydrodynamics are thus reasonable assumptions for building simple analytic 
models on how the IGM is reionized. If we further simplify by assuming an isolated point source 
in a homogeneous IGM and a constant temperature in the ionized region we get the well known 
Stromgren sphere model (Stromgren 1939). In this case the radius of the spherical ionized region 

1 /3 " 3/V " ^ 

grows as: r(t) = r s [I — exp(f/z rec )] where r s = 47ta ^ n 2 i s the asymptotic radius of the 
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ionized sphere and N-, is the number of ionizing photons produced per second by the source. From 
this model we can estimate that the typical velocity of ionization fronts in the IGM during the EoR 
is > 10 3 km.s” 1 . Let us emphasize that during the EoR, the luminosity of proto-galaxies typically 
increases on time scales much shorter that t rec , and thus the asymptotic Stromgren radius is never 
reached. 

Generalized solutions for the Stromgren sphere in an expanding universe have been formulated 
(Shapiro & Giroux 1987). In view of the future attempts to detect individual bubbles in the SKA 
tomographic data, it is also important to notice that the apparent shape of the bubbles along the 
light cone is not isotropic (Yu 2005; Majumdar et al. 2011). 

3. Ionization processes in primordial galaxies 

3.1 Escape of the ionizing continuum 

The ionization history of the IGM is regulated by the number of ionizing photons escaping 
from primordial galaxies. From an observational point of view this number can be inferred from 
the galaxy luminosity function at various redshifts. The most recent data give us indications to 
the amplitude and shape of the luminosity function up to z = 8 (e.g. Bouwens et al. 2014). But 
then modeling comes in to connect the luminosity function to the mass function, through a M/L 
ratio that may depend on a number of parameters. For an in-depth understanding of the physics of 
primordial galaxies however, it is useful to break the M/L parameter in two parts: the star formation 
rate (SFR) and the radiation escape fraction / esc , and more specifically for our purpose the escape 
fraction of the ionizing continuum. Indeed, primordial galaxies host dense gas clouds with short 
recombination time that will consume a large fraction of the emitted photons before they ever reach 
the slowly-recombining diffuse IGM. Dust may also play a role in absorbing photons during the 
later phase of the EoR. In any case, / esc and its dependence on the host galaxy and its environment 
is the result of complex radiative processes. 

From an observational point of view, measuring / esc during the EoR presents a challenge since 
the photons that do escape galaxies will be absorbed in the IGM and participate in the global 
reionization process, and thus will never reach us. Indirect probes can be used, such as the de¬ 
termination of the average photoionization rate in the IGM from Ly-a forest observations, or the 
Lyman-a emitters (LAE) luminosity since the intensity of the Lyman-a line is proportional to the 
recombination rate in the galaxy and thus to the amount of ionizing radiation that does not escape. 
Using such observations, a substantial level of modeling is needed to yield / esc . Observations pro¬ 
vide some constraints for the relative escape fraction between ionizing and non-ionizing photons at 
z ~ 3 (e.g. Steidel et al. 2001; Iwata et al. 2009). A few studies focused on higher redshift provide 
only tentative results: an upper limit of 0.6 at z=5.6 (Ono et al. 2010) or an inferred value of 0.3 
at z=6-8 (Finkelstein et al. 2012). However, there is a general agreement based on models that, 
if stellar sources are to be the main contributor to reionization, the escape fraction must increase 
toward higher redshifts and/or be higher in low mass galaxies (e.g. Alvarez et al. 2012; Haardt & 
Madau 2012; Kuhlen & Faucher-Giguere 2012; Mitra et al. 2012). 

Another way to produce estimates of / esc is to run numerical simulation that preferably include 
the coupling between the ionizing continuum radiative transfer and the dynamics. These simula¬ 
tions need to resolve to some extent the interstellar medium while preserving the cosmological 
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environment and have become possible only recently (Gnedin et al. 2008; Wise & Cen 2009; Ra- 
zoumov & Sonmier-Larsen 2010; Yajima et al. 2011, 2014; Paardekooper et al. 2013; Wise et al. 
2014; Kinim & Cen 2014). They do not use similar resolution nor similar modeling, and thus do 
not quantitatively converge in terms of the mean value and scatter of / esc as a function of redshift 
and host halo mass. One common trend seems the emerge however: / esc decreases with the mass 
of the host halo (with Gnedin et al. (2008) finding the opposite trend however, and Kinim & Cen 
(2014) finding almost no dependence), and increases with redshift (weak opposite tend in Kinim 
& Cen (2014)). This seems to indicate that low mass galaxies made the largest contribution to the 
production of ionizing photons, at least at high redshift, before radiative feedback quenched their 
star formation rate. 

3.2 Feedback on primordial galaxies 

The second important aspect of ionizing radiative processes in primordial galaxies is the feed¬ 
back on the gas (and thus star formation). This feedback is enacted through the photo-heating and 
photo-evaporation of the gas in galaxies and of the accreting gas clumps and filaments, but also to 
some extent through radiation pressure. 

In a forming galaxy made of near-primordial gas, cooling processes below 10 4 K are rather 
inefficient (low metal content, uncertain PF formation). Thus photo-heating the cold gas back to 
~ 10 4 K through ionization may substantially decrease the SFR. This is easily understood in < 10 8 
M 0 galaxies that have a virial temperature below 10 4 K: the gas will be expelled from the galaxy 
(Petkova & Springel 2011; Dale et al. 2012; Hasegawa & Semelin 2013). But it may also be true in 
more massive galaxies where the gas is not expelled but merely prevented from fragmenting. Such 
an effect is found by Hasegawa & Semelin (2013) for example. Such small scale processes can be 
implemented in large scale reionization simulations using sub-grid recipes (Iliev et al. 2007) 

On top of the feedback mediated by thermal pressure, radiation pressure may play a role in 
regulating star formation in primordial galaxies. Wise et al. (2012) observed, using numerical 
simulation, a significant quenching of the star formation rate when ionizing radiation pressure 
is included. This result has not been reproduced yet by other teams using different numerical 
techniques. It is definitely worth further investigation. 

Last but not least, the feedback from supernovae explosions definitely regulates the escape 
fraction. Although not a pure radiative process, it shapes the gas distribution around the sources 
and influences the escape fraction of ionizing photons (Kimrn & Cen 2014). 

The relative importance of these different processes probably depends on the mass of the host 
halo and is not the object of a consensus. Further numerical investigation is needed. 

4. Ionization processes in the IGM 

As we have seen in section 2.1, the mean free path of ionizing photons in typical IGM condi¬ 
tions during the EoR is very small at the ionization energy threshold compared to typical cosmo¬ 
logical distances. Consequently the ionization process sourced by UV photons takes the form of 
ionizing fronts sweeping through the IGM. The isothermal sound speed in the ionized gas behind 
the front is ~ 10 km.s -1 , and much lower in the neutral gas ahead of the front. The ionization front 
velocity (set by the local ionizing flux and baryon density alone if recombination is neglected) is 
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typically larger than 1000 km.s 1 in the IGM at z ~ 10 (see section 2.3). Consequently the gas is 
unable to respond dynamically, and the front is an extreme weak-R type (R for rarefied, see Spitzer 
1978), with negligible compression of the gas behind the front. This is the main justification for 
not coupling radiative transfer to hydrodynamics in many EoR simulation. It is worth mentioning 
however, that the ionization front not only slows down but even stalls in structures with moder¬ 
ate overdensity (~ 30) due to recombinations (e. g. Sobacchi & Mesinger (2014) and references 
therein). 

While there is little doubt that stars contributed a significant faction of the ionizing photons 
through their UV continuum, we have few constraints yet on the contribution of X-ray sources 
(X-ray binaries, Supernovas, AGN and more) to the reionization process. Because of their long 
mean free path, a reionization powered by X-rays only would show a diffuse and fluctuating ion¬ 
ization field without any fronts. It is however unlikely that X-rays contribute as much as the UV 
continuum to the ionization budget. Consequently the most likely picture is that of an IGM with 
fully ionized regions around the sources, separated by a sharp ionization front from near-neutral 
regions, ionized at most a few percent by X-rays. An important effect of X-ray penetrating into 
the neutral IGM however is the heating by secondary electrons (Furlanetto et al. 2006; Shull & van 
Steenberg 1985). Indeed when an X-ray photon ionizes an hydrogen atom, the energy in excess of 
the ionization threshold goes into kinetic energy for the freed electron. This electron then collides 
with neighboring atoms, producing secondary ionizations, excitations and heating. Heating up the 
cold IGM by a few K only has a strong impact on the 21 cm signal. 

The presence of helium in the primordial gas has a moderate impact on the reionization of 
hydrogen. With first and second ionization energies at 24.6 eV and 54.4 eV respectively. Helium 
filters and hardens the spectrum of Pop III and Pop II stars. In practice the emissivity above 54.4 
eV is small enough that the reionization of Hell is negligible at z > 6. On the other hand, Hel is 
reionized mostly at the same time as HI, moderately slowing the pace of the process by consuming 
a small paid of the radiations (e.g. Ciardi et al. 2012). Behind the ionization fronts, the temperature 
of the ionized IGM is somewhat higher than if Helium is not considered. Overall, the impact of 
helium reionization on the 21 cm signal is small. 

5. Recombination processes in the IGM 
5.1 Minihalos 

The very first stars formed in cosmological minihalos, much smaller than present-day galaxies, 
with total masses below ~ 10 8 M Q . The efficiency of this process is however uncertain. Indeed, the 
Jeans mass criterion for star formation requires that the gas be extremely cold, with temperatures 
of just a few hundred K. At primordial composition, consisting almost exclusively of Hydrogen 
and Helium, the gas in such small halos can only cool through Hn molecular line cooling. These 
molecules are fragile, however, and easily dissociated by Lyman-Werner radiation. Moreover, 
minihalos arc photoevaporated by ionizing radiation (Shapiro et al. 2004; Iliev et al. 2005b), which 
is a complex process that results in screening of the intergalactic medium by them and slowing 
down the large-scale ionization fronts (Iliev et al. 2005a; Ciardi et al. 2006). 

Furthermore, due to their shallow gravitational potential wells, the minihalo abundance is 
influenced by heating of the IGM by e.g. X-rays, as well as modulated locally by baryon-dark 
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Figure 1: Maps of the ionization fraction color coded form 0. (blue) to 1. (red) in simulations modeling the 
effects of minihalos (left) or not (right). From Ahn et al. (2012). 


matter velocity offsets (Tseliakhovich & Hirata 2010). The relative importance of these effects is 
still highly uncertain and requires much better knowledge than we currently possess of the early 
sources of X-rays and soft UV radiation: their nature, abundance and clustering. Nonetheless 
different aspects of this physics have been modeled in both simulations (Ahn et al. 2009; Baek 
et al. 2010; Ahn et al. 2012) and semi-analytical work (Santos et al. 2010). Both X-rays and 
Lyman-Werner radiations have very long mean free paths, of order 100 Mpc or more, and the 
baryon-dark matter velocity modulation acts on similarly large scales, all of which could result 
in significant fluctuations in the local abundance of minihalos and their star formation rates (e.g. 
Fialkov & Loeb 2013). The first star's forming in cosmological minihalos likely dominated the early 
reionization, but eventually the Lyman-Werner radiation from those same stars suppresses further 
star formation, stalling the reionization history until the larger-mass, atomically-cooling halos take 
over (Ahn et al. 2012). This results in significantly different early reionization history compared to 
the case without nrinihalos, much higher electron-scattering optical depth for CMB photons and a 
different reionization geometry (see Figure [I]). 

5.2 Lyman Limit systems 

Lyman Limit Systems (LLS) are structures of any mass that are optically thick to the UV 
background (in practice with column density larger than 10 17 cm -3 ). They act as screens for UV 
photons traveling in ionized regions. Analysis of QSO spectra show that LLS limit the mean free 
path of UV photons to ~ 50 cMpc at z = 6 (Songaila & Cowie 2010). Estimates at higher redshifts 
are uncertain. The main consequence of LLS would then be to slow down the end of reionization 
(in terms of the average ionization fraction reaching 1), by limiting the effective number of sources 
that contribute to the local ionizing flux. 

LLS are difficult to take into account in full numerical simulations because it typically in¬ 
volves resolving the internal structure and internal ionization processes in dwarf galaxies. Subgrid 
models can be devised and implemented. Using semi-numerical simulations Crociani et al. (2011) 
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and Sobacchi & Mesinger (2014) find that LLS have an impact on the history and geometry of 
reionization and thus on the 21 cm signal power spectrum. Alvarez & Abel (2012) also implement 
the effect of LSS in simulations. 

6. Conclusion 

At first glance, the process of reionization is simple: ionized bubbles around primordial galax¬ 
ies with sharp, supersonic ionization fronts expanding into a neutral IGM mildly heated and weakly 
ionized by X-rays. It is however regulated by two complex processes: the production and escape 
of ionizing photons from primordial galaxies, and the presence of dense, small scale structures in 
the IGM that act as photon sinks (minihalos and LLS). Both of these processes have an impact on 
the pace and geometry of reionization. On the one hand observations provide limited constraints 
on these processes at z < 6 whose extrapolation during the Eopch of Reionization is risky. On the 
other hand, simulations do not offer a consensus on these processes, mainly because they involve 
small scales that are difficult to resolve while keeping a statistically significant volume. Either sig¬ 
nificant progress will be made in the next few years or the interpretation of the observations of the 
21 cm signal with the SKA will have to deal with models involving a large parameter space. 
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